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A study is presented of a two-dimensional frustrated and 
dimerized quantum spin-system which models the effect of 
inter-chain coupling in a spin-Peierls compound. Employ- 
ing a bond-boson method to account for quantum disorder in 
the ground state the elementary excitations are evaluated in 
terms of gapful triplet modes. Results for the ground state 
energy and the spin gap are discussed. The triplet disper- 
sion is found to be in excellent agreement with inelastic neu- 
tron scattering data in the dimerized phase of the spin-Peierls 
compound CuGeC>3. Moreover, consistent with these neutron 
scattering experiments, the low-temperature dynamic struc- 
ture factor exhibits a high-energy continuum split off from 
the elementary triplet mode. 

PACS: 75.10.Jm, 75.40.Gb, 78.70.Nx 



I. INTRODUCTION 



The recent discovery of the anorganic spin-Peierls 
compounds CuGeC>3 jj] and a'-NaA^Os [||,|| has great- 
ly stimulated interest in low— dimensional magnetism. 
While many properties of these novel materials can be 
described in terms of quasi one-dimensional (ID) quan- 
tum spin-systems, clear evidence for a substantial de- 
gree of iwo-dimensionality (2D) of their magnetism has 
been found - most noteworthy by inelastic neutron scat- 
tering (INS) which displays a sizeable transverse disper- 
sion of the magnetic excitations in CuGeC>3 • In the 
present study I will establish a simple framework to in- 
terpret the spin dynamics of a frustrated and dimerized 
2D quantum spin-model with a particular focus on the 
low-temperature phase of CuGeC>3. 

CuGeOa is an anorganic spin-Peierls system with a 
lattice dimerization transition at a temperature T$p — 
14K Its structure comprises of weakly coupled 

Cu02 chains along the c-axis, with copper in a spin- 
1/2 state ||[l0). The nearest-neighbor (n.n.) exchange- 
coupling between copper spins along the Cu02 chains 
is strongly reduced by almost orthogonal intermediate 
oxygen states |rjj| . Therefore, next-nearest-neighbor 
(n.n.n.) exchange in CuGe03 is relevant. Both, n.n. and 
n.n.n. exchange, are antiferromagnetic (AFM) flTl| , [l2f im- 
plying intra-chain frustration. In addition, n.n., as well 
as n.n.n. inter— chain exchange is present which proceeds 
via the 02 sites |ll] , |l3| ] . This exchange is believed to be 
one order of magnitude less than the intra-chain cou- 
pling [Q^,[ll],[l3| and comparable to T$p. Therefore the 
inter-chain coupling should be relevant in the dimerized 
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FIG. 1. The J-X-a-fl-/3 model. Line segments refer to ex- 
change couplings for spins located at segment vertices. In- 
ter-chain couplings are shown for a single dimer site only, 
b = bet and c = ce c are the primitive vectors. 



phase. This may be a key element in the INS ^-01 an d 
magnetic Raman scattering Q |l6| data. 

A minimal model of CuGe03 which includes intra-, as 
well as inter-chain interactions is the J-\-a-u-{3 model 
0,|3 depicted in fig. |. The various line segments label 
the coupling strengths J, JA, JXa, Ja, and J/J./3 between 
the spins located at the vertices in this figure. J refers to 
the 'strongest' or dimer-bond the left vertices of which 
form the dimer lattice 1 € T>. Most important the dimer- 
ization in fig. [l] is staggered along the b-axis. This is 
realized both, in CuGe03 [[l0| as well as in o/-NaV205 
ifLsf , and turns out to be relevant for the magnon dis- 
persion. In fig. [j] an additional, so-called 'natural' label- 
ing of the intra-chain parameters is introduced, i.e J\, a, 
and 8. This notation is frequently used in the context 
of the ID dimerized and frustrated spin-chain limit. In 
CuGeC>3 Ji is approximately \QQK [ ^9pC| ]. Consensus 
on the precise magnitude of the intra-chain frustration- 
ratio a is still lacking. Studies of the magnetic suscepti- 
bility, which has been compared only to ID models, have 
resulted in a « 0.24 JT|] as well as in a « 0.35 @. 
This would place CuGeOs in the vicinity of the critical 
value a c ~ 0.2411 for the opening of a spin gap solely due 
to frustration ^]J. 6 resembles the lattice dimerization 
which is finite for T < T$p only. Values for the zero- 
temperature dimerization S(T = ) rang ing from 0.21 to 
0.012 have been suggested pi| , |T9| , po| , p2[ . Knowledge on 
the magnitude of fj, and u(3 is limited to |/x|,|/U/3| -C 1 

0H- 

Magnetic excitations in CuGeOs are clearly distinct 
among the uniform (U), i.e. T > T$p, and the dimerized 
(D), i.e. T < Tsp, phase. While the dynamic structure 
factor exhibits a gapless, c-axis dispersive two-spinon 
continuum similar to that of the ID Heisenberg chain 
above T$p p3| , well defined magnon-like excitations with 



1 



sizeable c- and b-axis dispersion have been observed be- 
low Tsp 0-^,0. These magnons are gapful and are 
split off from a continuum which, at zone-center, starts 
at roughly twice the magnon gap |p4-26,17||. 

The aim of this work is to study the magnetic proper- 
ties of the J-X-a-/j,-/3 model with a focus on the D-phase 
of CuGeC>3. First I will describe a bond-spin represen- 
tation of the J-X-a-/j,-f3 model which, in turn, is treated 
by an appropriate linearization. Next, results for the 
ground state energy, the spin gap, the magnon disper- 
sion, and the dynamic structure factor are contrasted 
against other theoretical approaches and are compared 
with experimental findings. Finally, details of an alter- 
native mean-field approach using the bond-spin repre- 
sentation are provided in appendix [X[ 



II. BOND-OPERATOR THEORY 

In this section the properties of the J-X-a-fj,-/3 model 
are discussed by representing the site-spin algebra in 
terms of frond-spin operators First, the essential 

features of these operators are briefly restated. Consider 
any two spin-1/2 operators Si and S2. The eigenstates 
of the related total spin are a singlet \s) and three triplets 
|f Q ) with a = x, y, z. These can be created out of a vac- 
uum |0) by applying the bosonic operators s' and v a 
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where [s,s^] = 1, [s (t) ,4 t} ] = 0, and [t a ,t\] = 6 a p. The 
action of Si and S2 on this space leads to the represen- 
tation 
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for the individual spin operators. Here £ Q/ 3 7 is the Levi- 
Civita symbol and a summation over repeated indices 
is implied hereafter. The upper(lower) subscript on the 
lhs. of (0) refer to upper(lower) sign on the rhs.. The 
bosonic Hilbert space has to be restricted to the physical 
Hilbert space, i.e. either one singlet or one triplet, by the 
constraint 
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Using (g) and ||) it is simple to check, that Si and S 2 
satisfy a spin algebra indeed and moreover that 
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In order to transform the J-\-a-fi-f3 model into the bo- 
son representation a particular distribution T of bonds, 
i.e. pairs of spins Sn and Si 2 , has to be selected. Here 



this selection will be based on the limit of strong dimer- 
ization (A, Act*) — * (0,0), or equivalently (a, 8) — > (0,1), 
and small inter-chain coupling (//, — > (0,0). In this 
limit the ground state is a product of singlets on each 
dimer-bond while the elementary excitations arc com- 
posed of the corresponding localized triplets. Therefore 
it is natural to place the singlet and triplet bosons onto 
the dimer bonds, i.e. Sn = Si and Si 2 = Si +C with 
1 £ T>. The transformed Hamiltonian reads 
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where each local Hilbert space is subject to the constraint 
(||) and, if not explicitly stated otherwise, the unit of 
energy is J hereafter. The inter-dimer matrix elements 
can be obtained from fig. [j] 
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where li,2,3 are defined in fig. [I] and t\ = A(l — 2a), 
t 2 = m(1 - 2/3), h = A(l + 2a), and U = /x(l + 2/3). 
As anticipated, the inter-dimer matrix elements a(l, m), 
6(1, m), and c(l, m) vanish in the strong dimer limit leav- 
ing the Hamiltonian diagonal in si and t\ a . 

In order to treat the local constraint and the dimer in- 
teractions approximations have to be made. To this end I 
will employ the Holstein-Primakoff (HP) representation 
of the bond-operators which has been detailed in [||| [3(j| . 
In this representation the constraint is treated by elimi- 
nating the singlet operator via s\ = s l = (1 — £icA Q )~ 1/ ' 2 - 
Moreover, after inserting this into the Hamiltonian only 
terms up to second order in in the triplet operators are 
retained. The latter procedure is analogous to the linear 
spin-wave approximation in systems with broken spin- 
rotational invariance. The linearized (LHP) Hamiltonian 
is given by 
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ek = — - [ti cos(2fc c ) + 2t 2 cos(fe&) cos(k c )] 
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FIG. 2. Ground state energies in the dimerized-chain 
limit: LHP-theory (solid) versus exact diagonalization |3l]] 
(solid dots), errors are less than marker size. Inset: small S 
limit. E*{5)/h = -[(l + 5)E g ((l-S)/(l + S),0)- Ef ethe ]. 
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FIG. 3. Ground state energies in the frustrated-chain 
limit: LHP-theory (solid) versus exact diagonalization ^] 
(solid dots), and DMRG (crosses), errors are less than 
marker size. E*(a)/Ji = -[E g ((l - 2a), 0)- E* ethe ]. 



where D is the number of dimers and k is a momentum 
vector restricted to a Wigner-Seitz cell B of the reciprocal 
lattice. For comparison with experimental data B is ori- 
ented with respect to the non-dimerized system, instead 
of the Brillouin zone of the dimer lattice, i.e. k = k c ) 
with bki, = 0...2-7T and ck c — 0...7T with b,c set to unity 
hereafter. is a a spinor with ^ Q = [ij^ £_ kc J and 
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Eqn. (ffl) describes a threefold degenerate set of disper- 
sive triplets. The triplets are renormalized by ground- 
state quantum-fluctuations. These are produced by the 
terms of type t^ a t_^ a and their hermitian conjugate. 
The excitation spectrum E^ follows from a Bogoliubov 
transformation 



Hlhp = -\d+ J2 ^(4 Aa + ^) (9) 

keB.a 

£ k = VI + 2e k , (10) 

where a k ^ are the Bogoliubov quasi-particles which are 
given by 
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where is a a spinor with $ kQ , = [i 

To conclude this section, I note that instead of treating 
the constraint by means of the HP representation one 
may also apply the so called bond-operator mean-field 



theory (MFT) of [§7| . The application of this method to 
the J-X-a-fi-(3 model is d etailed in appendix |A|. As will 
become evident in section [II B this technique seems less 



well suited in the present context. 



III. RESULTS 

In the following sections the consequences of the LHP 
representation of the 2D J-X-a-fi-(3 model will be con- 
trasted against other known results in the limiting case 
of the ID Ji-a-5 model as well as INS data observed on 
CuGe0 3 . 



A. Ground State Energy 

From (fj) it is obvious, that the ground state energy per 
lattice site, E g , is equal to —3/8 at the point of complete 
dimerization. This is the proper energy gain for a bare 
singlet formation. For arbitrary t\ and t% 

E g {hM) = -| + A I dk c {[l-2t 2 cos(k c ) (12) 
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where E is the complete elliptic integral of the second 
kind. For vanishing inter-chain coupling this simplifies 
to 



E g (t u 0) = -9/8 + — VT^hEi 
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which can be compared to existing results in various re- 
gions of the (<5, <5)-plane of the frustrated and dimerized 
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FIG. 4. Spin gaps in the dimerized-chain limit: LHP — 
theory (solid) versus 3rd-order perturbation theory ( 
dashed), MFT (dotted), and exact diagonalization ]3l| (solid 
dots), errors are less than marker size. 




FIG. 5. Spin gaps in the frustrated-chain limit: LHP — 
theory ( solid ) versus 3rd-order perturbation theory 17] 
(dashed), MFT (dotted), and two DMRG calculations j 3E ] 
(crosses), and [M (solid dots), errors are less than marker 



spin-1/2 chain. In particular at the isotropic Heisenberg 
point, i.e. (5, 5) = (0,0), E g = 3/ (V2n) - 9/8 « -0.4498 
which agrees reasonably well with the Bethe Ansatz re- 
sult E^ ethe = l/4-ln(2) fa -0.4431. In fig. |the ground 
state energy along the (a = 0, <5)-line, i.e. for a dimerized 
chain, is contrasted against results from exact diagonal- 
ization J3lj . Deviations are within the numerical error of 
the diagonalization data for S > 0.05. For dimerizations 
below 0.05 small differences are caused by an unphysical 
extremum in the LHP energy which is visible in the inset. 
While the agreement is encouraging along the (a — 0, 6)- 
line, only qualitative consistency is to be expected along 
the (a, 6 = 0)-line since this region is more distant from 
the strong dimcr limit (a, 5) ~ (0, 1). This is shown in 
fig. which compares E g (a) with exact diagonalization 
[32 1 and density-matrix renormalization group (DMRG) 
[33 1 data for the case of a frustrated chain. This figure 
demonstrates that the LHP approach overestimates the 
frustration induced loss of ground state energy at inter- 
mediate a. 



B. Spin Gap 

The LHP approximation does not break the spin- 
rotational invariance and leaves the system in a quan- 
tum-disordered ground state. This is consistent with a 
spin gap A of the triplet dispersion which, for positive t\ 
and £2, is given by A = y/l — t± — 2t2 and is situated at 
k = (0, 0) and (it, it). If ti and ti are such that the triplet 
modes turn massless, i.e. A = 0, a quantum phase- 
transition towards antiferromagnetism (AFM) will occur 
at T = 0. In terms of a and S, the AFM instability-line 
is located at a + S — £2(1 + 8). Beyond this line the LHP 
approach breaks down. 



Next the LHP spin gap at ti = is compared to known 
results for the frustrated and dimerized spin-1/2 chain. 
In fig. [I| A (5) is contrasted against findings of exact diag- 
onalization |3l| ] , perturbation theory up to third order in 
A and a solution of the MFT equations (|A5|-|A7|) 

which are discussed in appendix This figure displays 
reasonable agreement between the first three of these ap- 
proaches for all values of 5. In addition it shows that the 
MFT suffers from the inability to close the spin gap at 
the isotropic Heisenberg point |37||3^|. This caveat ren- 
ders the MFT unsuitable for the case of a systems with 
nearly massless spin excitations, e.g. CuGeC>3. 

As noted in the previous section, analytic approaches 
based on the strong dimer limit are less reliable if con- 
sidered along the (a, S = 0)-line. Nevertheless, a com- 
parison of A (a) as obtained from the LHP theory with 
various other techniques is instructive and is shown in 
fig. H Qualitatively, all analytic methods depicted ex- 
hibit a tendency of the spin gap to increase as a increases 
J35| however agreement with the DMRG data J|J|[36| is 
absent. In particular, while 3rd-order perturbation the- 
ory and MFT show no critical frustration-ratio a c and 
have a finite gap for all < a < 0.5 the LHP represen- 
tation leads to a critical frustration a c = 0. Quite re- 
markably there are also substantial differences between 
the two DMRG results for a > 0.5 Estl. 



C. Triplet Dispersion 

In this section the relevance of the model with respect 
to the spin excitations in CuGeOa will be assessed by 
comparison of with INS data for the D-phase. In 
particular the role of two-dimensionality will be consid- 
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ered. The essential effect of a finite hopping amplitude 
t-2 is a mixing of the b- and c-axis dispersion. This mix- 
ing is due to the staggering of the dimerization along 
the b-axis and leads to the ti cos(fc{,) cos(/c c )-term in et- 
Thus, i?k involves terms of different periodicity in k c , i.e. 
cos(2fc c ) and cos(fc c ). Therefore, in contrast to the quasi 
ID case, the degeneracy of the triplets at the momenta 
(0, 0) and (0, tt) is lifted. An identical reasoning based on 
the first-order contribution of a 3rd-order perturbation 
theory has been given in |l7| . Expanding E^ in terms 
of ti and t 2 I find agreement up to first order with the 
dispersion given in |17| ]. 

In fig. | INS data of the magnon dispersion in CuGeC>3 
are shown for momenta along the edges of the reciprocal 
lattice cell from k = (0, 0) to (0, tt) as well as from (0, 0) 
to (tt/2, 0) Although data exactly at (0, tt) is lacking, 
it seems very likely from this figure that the triplet exci- 
tations in CuGeC>3 are not degenerate at (0, 0) and (0, 7r). 
In order to demonstrate that the J-\-a-/i-{3 model can 
account for the observed dispersion fig. || contains a com- 
parison of i?k with the INS data. Rather than performing 
this comparison by a least-square fit, has been identi- 
fied with the magnon energies only at k = (0, 0), (0, 7r/2), 
and (7r/2,0). This fixes J = 11.5meV, h = 0.859, and 
t 2 = 0.054 unambiguously and limits the b- to c-axis 
coupling-ratio to roughly 6% which is consistent with 
P-p|,^3|. While the preceding leaves a, S, f3, and /j, unde- 
termined, a can be fixed using an additional input, i.e. 
6 = 0.012. This is within the range of values suggested 
in the literature, i.e. 5 = 0.21. ..0.012 and 
implies natural parameters of J\ = HAmeV = 132i^ 
and a = 0.059. Allowing for a larger dimerization leads 
to smaller Ji and a. In particular a = with J\ = 12AK 
is reached for <5 = 0.076. Parameters with a slightly 
smaller (larger]intra(inter)-chain exchange have been es- 
tablished in (n), i.e. J x = 9.8meV = 114K, 5 = 0, and 
t 2 = 0.12 however with 5 = 0.12. 

The agreement displayed in fig. ^ is satisfying and 
demonstrates the main point of this section, i.e. that 
diagonal triplet-hopping can account for the observed 
asymmetry of the INS data in the D-phase of CuGc03. 
Moreover, the values of a obtained suggests that the 
intra-chain frustration in CuGeOs is significantly smaller 
than that derived from the purely ID Ji-a-5 model 
HQ, i.e. 0.24 < a < 0.36. The determination of 
model parameters however is in need of further studies to 
improve on the quantitative reliability of the frustration- 
dependence of presently available approaches for the 2D 
case. This is in contrast to fl7jj where quantitative evi- 
dence for a ~ 0, based on 3rd-order perturbation theory, 
has been suggested. 



D. Dynamic Structure Factor 

Magnetic excitations are observed by measuring the 
dynamic structure factor S^q, to) which is related by an- 




FIG. 6. INS data of c-axis (upper panel) and b-axis 
(lower panel) dispersion of the triplet mode in the D-phase of 
CuGe0 3 at T = 1.8K |J (solid diamonds) versus LHP-theory 
(solid) for J = ll.5m.eV, U = 0.859 and t 2 = 0.054. ( 
b = c=l). 



alytic continuation and the fluctuation dissipation theo- 
rem S(q, oS) = Im[x(q, w)]/(l — e~ w / T ) to the dynamic 
spin susceptibility 



X^(q ) r)=(T T [^(r)^]) 



(14) 



Here q is the momentum, r the imaginary time, and 
T T refers to time ordering. Since the physical spin is a 
composite operator of the bond-bosons the information 
obtained from the dynamic susceptibility is not restricted 
to the triplet dispersion - even at the LHP level. This 
will be clarified in the remainder of this section. 

Within the LHP representation the spin operator in 
momentum space is 



S$ = J(l-e^)(*U+*- q J 
1 
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Two qualitatively different excitations appear on the rhs.: 
a sharp triplet mode due to the first term and a contin- 
uum of two-triplet states due to the second. The mo- 
mentum dependent form factors (1 =F exp(ig c ) lead to a 
vanishing weight of the the triplet mode (continuum) at 
q c = (q c = tt). Since Hlhp does not conserve the bare 
triplet number, the continuum is divided in between two 
excitations of different nature, i.e. virtual excitations 
at energies -Ek+q + and real excitations at energies 
i?k+ q — Ek. While the former are due to ground-state 
quantum-fluctuations and occur at all temperatures, the 
latter result from excitations across the spin gap and are 
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FIG. 7. Energy bounds, f2c(q), for both triplet continua 
(dashed dotted) versus triplet-mode energy, -B(q), (solid) as 
well as combined spectral weight of continua, Wc(q), (dotted) 
at T = 0.05 versus spectral weight of triplet mode, Wr(q), 
(dashed) for momenta along qb — q c - (b = c = 1, J = 1 and 
ti,t a as in fig. p[) 



1.0 



1 0.5 
o 

II 

cr 

"J= 0.0 
01 1.0 



3 

-S0.5 



0.0 









1 



FIG. 8. Triplet continua for various momenta along 
qt — q c and for two temperatures T = 0.05 (a) and T — 0.15 
(b). Bottom-most curve corresponds to q = with 0.1 
incremental y-axis offset for each consecutive momentum. 
Solid diamonds: bounds of triplet continua. Solid square: 
triplet-mode energy. (6 = c = 1, J = 1 and ti,ts as in fig. 



present only at finite temperatures. Evaluating the dy- 
namic susceptibility by standard methods I obtain 

Xa^q, ^n) = 6 a p- { (cos(g c ) - 1W — ^ 



(cos(<? c ) + l) — 



1 + e k+q + £k + -Bk+q-Bk . 
-E"k+q-Ek 

n(£ k+q ) - n(E k ) 

lUJ n + -Ek+q — ^k 

^1 + Ck+q + £k — £'k+q£ , k)(-Ek+q + £?k) 
^k+q-Bk 

n(£ k+q ) + n(£ k ) + 1 



(iLOn) 2 - (£ k+q + £ k ) 2 



. (16) 



where io n = 2mrT is a Bose Matsubara-frequency. 

As anticipated the dynamical susceptibility (|l^) ex- 
hibits finite spectral intensity at E^, -Eic+q — -E k , and 
^k+q + -E-k- This is summarized in fig. [7] which displays 
the two continua with respect to the triplet mode along a 
particular momentum space direction. Parameters iden- 
tical to those of fig. || have been chosen. The triplet 
mode is situated in a gap between the low-energy con- 
tinuum due to thermal excitations and that at high en- 
ergies due to quantum fluctuations. At the zone center 
the latter is gaped by 2A. This is consistent with recent 
INS data from the D-phase of CuGe0 3 ||. In addi- 
tion to the magnon these experiments indicate an 'unex- 
pected' high-energy continuum which is separated from 
the magnon by an additional gap. The observed zone- 
center continuum-to-magnon gap-ratio is approximately 



2. Deviations from the latter value of 2 are possibly due 
to triplet-triplet interactions which are beyond the LHP 
approach. Additionally, fig. (?] shows the weight of the 
triplet mode, as well as that of the combined continua at 
T = 0.05 J. At T <C J the continuum weight is almost 
completely due to quantum fluctuations while, indepen- 
dent of temperature, the triplet weight is given by 



W T (q) 



7T 1 — COs(g c ) (jr,n) 7T 

8 ~E a = 4A 



(17) 



The leftmost expression refers to the maximum of the 
triplet weight. The corresponding momentum, i.e. (7r,7r), 
indicates the instability towards AMF as A — > 0. 

Figure ^ depicts the spectral intensity Imfxa^q, oj n — > 
—iuj + v)] of the continua as a function of frequency for 
two temperatures and for various momenta along a di- 
rection in reciprocal space identical to that of fig. [7]. The 
k-sums of (|l6|) have been performed numerically on a 
600x300 lattice. In order to obtain sufficient smooth- 
ing the frequency has been shifted off the real axis b 
v = 0.05. The parameters correspond to those of fig 
As a guide to the eye the position of the triplet mode 
is labeled by solid squares while solid diamonds in this 
figure label the exact locations of the spectral bounds 
of the continua. Although smeared due to the imaginary 
broadening, van-Hove singularities are clearly observable 
at the spectral bounds as well as other characteristic ener- 
gies within the continua. Evidently the quantum fluctua- 
tions exhibit largest weight at intermediate wave vectors 
while at higher temperatures additional weight appears 
at smaller momentum due to thermal excitations. Even 
though the intensity in fig. || decreases both, as q ap- 
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proaches (0,0) and (w,tt), only in the latter case this is 
due to the form factor in ( |l5| ) while in the former case 
this is a consequence of the conservation of the total spin. 



IV. CONCLUSION 

In summary I have studied static and dynamic prop- 
erties of a frustrated and dimerized 2D quantum spin- 
model using the bond-operator method. The ground 
state energy and the spin gap have been gauged against 
known results from the ID limiting cases of this model. 
Effects of dimerization are found to be described almost 
quantitatively while the influence of frustration is cap- 
tured qualitatively. The dynamic structure factor has 
been analyzed and displays two characteristic features, 
i.e. a well defined magnon excitation and a temper- 
ature dependent continuum. The magnon dispersion 
shows a characteristic lifting of degeneracies, different 
from purely one-dimensional models, and agrees very 
well with INS data on CuGeC>3. The low-temperature 
continuum exhibits a zone-center gap twice that of the 
magnon. This is also consistent with INS experiments on 
CuGe0 3 . 
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APPENDIX A: BOND-OPERATOR 
MEAN-FIELD THEORY (MFT) 

An alternative approach to the Hamiltonian (^) arises 
by introduction of a set of local Lagrange multipliers 771 
to enforce the constraint (0) 



Hmft - 
1 x 



D(-^--s 2 



2 Z^*k« 



— 77s 
s 2 e k 



5 ^ 
s 2 e k 



(A2) 



s 2 e k 



s 2 e k 



k a 



with nota tion s equivalent to (Q) . Analogous to the latter 
equation (A2) represents three dispersive triplet excita- 
tions, however, with a modified dispersion relation. After 
a Bogoliubov transformation one gets p2fl 



H 



MFT 



/ 3 3 2 

D( s 2 
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77 s 
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MFT 
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(A3) 



(A4) 



where d = 2s 2 /(l/4- 77). The b ka quasi-particles follow 
from expressions identical to (Llj) with e k , 1 + e k , and i? k 



replaced by s 2 e k , \ — 77 + s 2 e k , and Eyf FT , respectively. 
The Lagrange multiplier and the singlet amplitude have 
to be determined by solving the saddle-point equations 
(dH^px/dr)) = and (dHMFr/ds) = 0. This leads to 







2 3^1 



de k /2 



D^ B VT+d^ 



rfe k 



«^L&k 3 



(A5) 
(A6) 



These selfconsistency equations can be solved by combin- 
ing them into a single one for the variable d only 



d = 5 

(T=0) 



D ^ JT 
ke£ v 



kei3 



3_ 

D 



1 



(A7) 
(A8) 



where rj follows by insertion of d into ( |A6| ). This com- 
pletes the description of the bond-operator MFT. 



H = H 



lev 



mislay +t\ a t la -l) 



(Al) 



To treat this Hamiltonian one replaces the local con- 
straint by a global one, i.e. r/\ = 77, and introduces a 
mean-field (MF) decoupling of all quartic terms leading 
to an effective quadratic Hamiltonian p7| . This Hamil- 
tonian has an overall negative prefactor to the SjSj-term 
which implies Bose condensation of the singlets. There- 



fore s 



s is assumed. Moreover it can 



be shown that contributions from the triplic and quartic 
triplet-terms H2 and H3 to the MFT can be neglected 
The MF Hamiltonian reads 
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